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ABSTRACT 

We derive a formula for the anisotropic density distribution around haloes and voids 
on large scales. Our model assumes that the orientation of non-linear structures is 
strongly correlated with the Lagrangian shear, and gives a qualitative understanding 
of the recent detection of an enhanced clustering signal along the major axis of haloes. 
We also show that the measured amplitude is inconsistent with a model in which the 
alignment is produced by the initial inertia rather than shear tensor. 
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1 INTRODUCTION 

The clustering of matter at late times provides important 
constraints on cosmological models. Our understanding of 
the signal is best on large scales, wher e it can be de scribed 
by perturbation theory well (see e.g. Peebles 1980). E.g., 
the Baryonic Acoustic Oscillations in the power spectrum 
(BAO), which appear as a spike in th e two-point corr ela- 
tion function, lie in this regime (see e.g. lEisensteinll2~005l ) . In 
addition to the simple correlation function, there are other 
ways to extract information from the matter distribution. In 
redshift-space distorted measurement s, the two-po int corre- 
lation funct ion is anisotropic (see e.g. Kaiser|[l987l or the re- 
cent work of lSchlagenhaufer et alj2012h . and this anisotropy 
can be used to constrain cosmological parameters. However, 
certain real-space measures of clustering are also expected 
to be anisotropic. E.g., galaxy clusters are typically triax- 
ial, and this triaxiality has long been know n to align with 
the s urrounding large scale structure (e.g., ISmargon et al.l 
|2012| . and references therein). On smaller mass scales, galaxy 
spins are also known to align wit h the environment (e.g., 
IJones et all bold : iLee fc PenllioOll ) . 

Recently, [Faltcnbac her et "all (|2012j ) showed that, in 
their numerical simulations of hierarchical gravitational 
clustering, the cross-correlation function between haloes and 
the surrounding mass was anisotropic: this correlation be- 
tween halo shapes and large scale structure extended even 
to the large scales relevant to BAO studies, and affected the 
zero-crossing of the correlation function. This motivates our 
work, which attempts to model this anisotropy. 

Our model, which is described in Se ction 2, is based 
on the assumpt ion that halo s hapes (ILee et alj 120051 ; 
iRossi et af1l201ll ) and orientations (|Lee fc Penll200ll ) at late 
times are correlated with the properties of the initial La- 
grangian field from which they formed. This is a fundamen- 
tal ingredient in models where haloes form from a triaxial 



collapse l|Bond fc Mversl Il996l ; ISheth et all l200lh . In such 
models, the Lagrangian deformation or shear tensor plays 
a key role, as its eigenvalues can be used to distinguish be- 
tween haloes, filaments, walls and voids. We illustrate our 
model for the two extreme cases: haloes and voids. A final 
section discusses potential applications and extensions of our 
work. 

Although our analysis is general, when we illustrate our 
results, we will assume a ACDM model with h = 0.73, 
fi cdm = 0.215, n bar = 0.035, Ov = 0.75, and a 8 = 0.9. 
These values facilitate co mparison with the simulations of 
iFaltenbacher et al.1 i|2012f ). 



2 ANISOTROPIC DENSITY DISTRIBUTION 
AROUND HALOES 

Despite the fact that haloes are highly non-linear ob- 
jects, their formation encodes information about the ini- 
tial (Lagrangian) fields from which th ey formed (e.g., 
iPress fc Schechter|[l97i : ISheth et al.ll200ll ). so, for example, 
one expects the shape, spatial orientation and sp in of a halo 
to be correlated with t he initial tidal field (e.g.. ILee fc Pen! 
120011 : IRossi et al.|[201ll ) . 



2.1 The shear 

The initial tidal field or shear tensor is defined as 



1 <9 2 $(q) 
cro dqidqj 



where A$(q) = <5(q) 



(1) 



where q is the Lagrangian spatial coordinate, $ is the La- 
grangian potential at q, and Oq ls the variance of the La- 
grangian density 8 smoothed on scale R. This variance de- 
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pends on the power spectrum P(k) and the smoothing filter: 
a) = / dkk 2 k 2i P{k)Wn{k). (2) 



For the ACDM parameters given earlier, we o btain the linear 
theory P{k) from CAMB (|Lewis et al.ll2000l ). The resulting 
dj decreases monotonically as R increases. As a result, in 
the excursion set description of haloes or voids l|Bond et al.l 
Il99ll ; [Sheth fc Tormenll2002l '). ao decreases as the halo mass 
or void radius increases. This will be important when we 
wish to relate our results to the halo-based measurements 
in simulations. 

In what follows, we would like to estimate the 
anisotropy in the correlation between the density distribu- 
tion at one position given that the shear at another position 
(which we will take to be the origin) satisfies some set of 
constraints. That is to say, we are interested in 

/ H d*p(0/d*(r)J(r)p(«(r)|£) 
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f s dtp(0 (5(r)|£) _ / s d£K6 (<5(r)|£) 
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P(H) 



(3) 



where £ = (£11, £ 22 , £33, £12, £23, £13) is a vector made from 
the components of the shear tensor and, for simplicity, we 
omit the distance argument r if the quantity in question is 
taken at the origin: e.g. £ = £(0). In the expression above, 
S is the region in £-space where the conditions on the shear 
field (associated with halo or void formation) are satisfied; 
we use P (H) to denote the integral over this region. For 
Gaussian initial conditions, p(£) is a multivariate Gaussian. 

2.2 Average of S(r) with conditions on the shear 

The Gaussianity of p(<5(r)j£) means that 

(5(r)|£} = (5(r)®£}(£®£)- 1 £ (4) 
£) nor (£ ® £) 1 depend on £. There- 



where neither (S(r 
fore 



(5(r)|E> = (<5(r)®£) (^^H) 



(5) 



The first two terms depend only on the correlations between 
S at one position and the shear tensor £ at another. Such cor - 
relations have been c omputed be fore Doroshkevichl (|l970T ) ; 
iBardeen et all i] 1986ft ; lDesiacauesl (|2008l ). We us ed the form 
for (£ g(r)£ki) that is given in the appendix of iDesiacquesI 
(2008). Although this expression can be worked out exactly 
for the 6x6 covariance matrix associated with £, it is sim- 
pler to work in the coordinate system which is aligned with 
the principal axes of the shear tensor. In this case, we find 
that 
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Similarly, using the fact that 5/ao = where ao was 

defined in equation {T}, we find 
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with f being a unit vector and 
1 



A„(r) 
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dkk 2 j n (rk)W R (k)P(k), 



(7) 



(8) 



where j n is a spherical Bessel function, and Wr is the win- 
dow function of scale R over which the shear is smoothed. 
Since 

a 2 A 2 (r) = £ 2pt (r)-£ 2pt (r) (9) 

where £ 2p t(r) = (Jq Ao(r) is the usual angle- averaged two- 
point correlation function, and 

ww=4 r dfr 2 £ 2pt (f) (10) 

r Jo 

is its volume average, equation ([7]) can be cast into a more 
intuitive form: 



(gtofa) 6 pt (r)-£ 2pt (r) , . Ur) 
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(11) 



In this expression, one should think of £ as the overdensity 
within r, and £ — £ as the difference between the overdensity 
within r and that at r itself. 

Inserting this in equation ([5]) yields 

(6(r,fi)\E) = (5\E)Ao(r) - 5 (£\S) A a (r) ( 12 ) 

where fi = cos 9, P 2 (/i) = (3/x 2 — l)/2 is a Legendre polyno- 
mial, 



(<5|S) = a (£11 +£22 + £33|H), 

/ — It £22+£33,-\ 

{t\a) = a (£u 1=) 



(13) 
(14) 



In this form, it is clear that the first term on the rhs 
of equation (|12[) . (S\S) £ 2p t(r)/cr 2 , is the spherical average 
of the full expression. Therefore, the prefactor should be 
thought of as a 'linear bias factor' 
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6(H) = <<5|H)/a 2 



(15) 



coming from the constraints S. (We provide an explicit ex- 
ample of this in the next section.) The angular dependence 
comes from the second term, which vanishes when £ = 0, 
i.e., for spherical objects. 

Since both 6(H) £ 2p t(»") and its angle average 6(H) £ 2pt (?") 
can be measured (indeed, these are the traditional 2-point 
measurements), our model can be written as 



(S\E) A (r) 
where 



A$S) = b(l\E)/(6\3). 



(16) 



(17) 



The left hand side is the ratio of observables, and the right 
hand side shows that it is the product of a scale- independent 
amplitude, and a separable function of scale r and angle. 
Since the r dependence is completely specified by measur- 
able angle- averaged quantities, and the /j, dependence is sim- 
ply that of a quadrupole, the amplitude is the only free pa- 
rameter in our model. In this respect, the expression above 
should be thought of as providing a generic fitting formula 
with just one free parameter, the amplitude. In our model, 
this amplitude encodes information about the alignment be- 
tween the tracer field and the large scale environment. 

2.3 Illustrative constraints 

The averages over H in equations (|13[) and (|14[l can be calcu- 
lated for many scenarios. To gain intuition, suppose that we 
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Figure 1. Equation i l I with J^i&i > 1.66/ao an d §22,33 > §11 > 0.41/cto for different mass scales. 



identify haloes with regions in the initial field for which a ll 
three eigenvalues were positive fe.g. iLee fc Sh andarin 19981 ). 
In this case, requiring 



S = {£33 > £22 > £ n > 0} 



means that 
(fu|S> = 
{62 18> = 

(63 1 3) = 
where 
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(19) 
(20) 
(21) 

(22) 
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and A(H) « —7/4. This shows that the bias factor increases 
as <7o decreases. We remarked earlier that, in the excursion 
set approach, large masses have small ao. Therefore, this 
model has the monopole part of the signal increasing as mass 
increases, but the ratio of the monopole to the quadrupole 
is independent of halo mass. 



2.4 More realistic conditions for haloes and voids 

Figure [1] shows the result of evaluating equation {12} numer- 
ically with S given by the requirement that J]V > 1.66/cto 
and £22. 33 > £u > 0.41/cto for a range of choices of ao. These 
constraints on the £y were motivated b y the spher i cal co l- 
lapse model and additional analysis in lLam et all (2009). 
The value s of op were chosen to match the halo masses 
quoted bv lFaltenbacher et al.l (|2012l ) in their analysis of the 
anisotropic clustering around haloes in simulations. 

It is encouraging that the amplitude and scale depen- 
dence of our predicted signal are qualitatively rather similar 
to theirs, although our angular dependence is slightly larger 
than they report. While errors in measuring the direction of 
the long axis of each halo will weaken the measured signal, 
thus alleviating some of the discrepancy, there are at least 



three other reasons why this may happen. First, our model 
of the Lagrangian patches which become haloes is crude, and 
may be inadequate. Second, h aloes are not, in f act, perfectly 
alighted with the shear (e.g. ILee fe Per]|200ll ). And third, 
our model is based on linear theory; recent work has shown 
that nonlinear evolution will induce a quadrupole even if 
none is i nitially present, a nd will modify it if there is one 
initially l|Chan et al.ll2012l . although, on the scales of most 
interest here, this is expected to be subdominant). Account- 
ing for this will weaken the predicted signal. 

Other non-linear structures include filaments, sheets, 
and voids. In triaxial models, this classification is related to 
the eigenvalues of the shear (|Shen et al.ll2006T h All positive 
eigenvalues describe a halo, one negative gives a filament, 
two negatives a sheet, and all three negative a void (e.g. 
lHahn et al.ll2007l ). For voids a reasonable set of conditions 
on the eigenvalues is > £22,33 > £11 and J2i < — 2.8/ao. 
Again, the latte r condition is from the spherica l evolution 
model (see e.g. ISheth fc van de Weygaertl l2004f i. Figure [2] 
shows the result. In the direction of the major axis, there is 
a positive boost similarly to the case of haloes. Of course, 
since voids are driven towards spherical symmetry, their ori- 
entation may be harder to measure accurately. 



2.5 The shear vs the inertia tensor 

Haloes have been shown to align with the Lagrangian shear 
l)Dubinskill992l ; rLee fc Penll200lh . and tri axial collapse mod- 
els relate halo shapes to the initi al shear (ISheth et al.ll200ll ; 
ILee et allfeoOSl i lRossi et al.ll201lh . These findings single out 
the shear tensor from others that could potentially define 
the axes of structures. An alternative to the shear is the ini- 
tial inertia tensor (the matrix of second der ivatives of the 
densi ty field). E.g., in the peaks model of ()Bardeen et al.l 
1986), haloes form at the peaks of the Lagrangian density 
field so their shape is given by the inertia tensor. 

The matter distribution around a peak with condi- 
ti ons on the inertia t ensor is given by equation (7.8) 
of iBardeen et all (|l986h . Angle averaging their expression 
around axis 1 yields 



- yx 



-5(Ai 



A (r) 
A 2 + A a 



7^ - 



■A„(r) 
)£a(r)ft(A0, 



(24) 



4 P. Papai & R. K. Sheth 




Figure 2. The same as Figure[JJbut with < — 2.8/o"o an d > §22, 33 > §11. 



where Ai, A2 and A3 are the eigenvalues of ~d 2 5(0)/dqidqj , 
X = X)i A,, V = 5/<7o, 7 = o"i/(oW2), and 

A„(r) = — ^ / dfcfc 4 j n (rfc)M/ fl (fc)P(fc). (25) 

27i- J o- O 0"2 J 

It is worth noting that Ai — (A2 + As)/2 is (3 y + z)/2, where 
y and 2 are the anisotropy parameters of iBardeen et al.l 
( 198^ 1. 

Although we have yet to average over the Ai, compar- 
ison with equation (|12|l shows that the first two terms in 
equation (124[) will yield the monopole, and the final term a 
quadrupole. Although the quadrupole here depends on the 
eigenvalues A of the inertia tensor in the same way that the 
quadrupole in equation (|f 2|) depends on the eigenvalues of 
the shear tensor, we might expect the amplitude here will 
be much smaller. This is because the integral which defines 
A2 has two additional powers of k compared to that which 
defines A2. However, we must also check that the average 
over the Ai does not yield a large amplitude to compensate 
for this difference. 

To see that this will not happen, note that on large 
scales, the leading order contribution to the monopole is 
given by the first term on the rhs of equation (|24[) . With 
the peak constraints, the quantity v = (y — 7x)/(I — 7 2 ) 
is a Gaussian variate, and it is independent of the Ai. 
This means that we can choose the constraints on v to 
match the monopole of equation l|12fl. leaving us to per- 
form an independent average over the dist ribution of pa- 
rame ters x, y, and 2 (e.g. Appendix C of IBardeen et al.l 
1986). The resulting angular dependence is so weak that 
it cannot be distinguished from the solid curve shown in 
Figure CD Since th e nonlinear evolution effects described by 
IChan etalJ (I20I2T) cannot, on their ow n, account for the 
signal seen by iFaltenbacher et al.l l|2012f ). we conclude that 
the initial shear tensor matters for the angular dependence, 
whereas the inertia tensor does not. 



3 DISCUSSION 

In this paper, we calculated the anisotropy in the linear den- 
sity field when conditions are placed on the Lagrangian shear 
field. If the shear field is strongly correlated with the shapes 
and orientations of nonlinear haloes, then this calculation 
should be closely related to the anisotropy of the halo-mass 



cross-correlation function, which is most easily seen when 
the mass field around haloes is stacked after aligning alo ng 
the major axis of the halo (e.g. IFaltenbacher et al.ll2012T ). 

For haloes, our model (equation If 2[) captures the main 
features of the cross-correlation function measurement. The 
signal along the long axis is stronger than perpendicular to 
it, but produces a less prominent BAO feature (Figure [TJ. 
We predict a similar effect for voids (Figure [2]). The overall 
amplitude of the effect is slightly stronger in our model than 
in the simulations. This may be due to inadequecies in our 
crude model which relates halo formation to the local shear; 
it may arise because the alignment between haloes and the 
shear is not perfect; or nonlinear evolution may have had 
a small effect (see Section I2,4|l . We also argued that the 
measured amplitude is inconsistent with a model in which 
the aligment is produced by the initial inertia rather than 
shear tensor (Section 12. 5p . 

Absent a model for halo or void formation, our equa- 
tion (|I2|I may be treated as a one-parameter family which, 
given the spherically averaged measurement describes the 
anisotropy. This parameter depends only on the local shear, 
and so may be used to constrain models of halo formation 
and alignment. Further work can be done to incorporate 
redshift distortions and nonlinearities into the model. Also, 
tests on simulations are needed in order to identify system- 
atics that can affect the validity of the model. Finally, we 
are in the process of checking if this sort of measurement 
can yield useful constraints on modified gravity models. 
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